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We study Doppler effects in curved space-time, i.e. the frequency shifts induced on elec- 
tromagnetic signals propagating in the gravitational field. In particular, we focus on the 
frequency shift due to the bending of light rays in weak gravitational fields. We consider, 
using the PPN formalism, the gravitational field of an axially symmetric distribution of 
mass. The zeroth order, i.e. the sphere, is studied then passing to the contribution of 
the quadrupole moment, and finally to the case of a rotating source. We give numerical 
estimates for situations of physical interest, and by a very preliminary analysis, we argue 
that analyzing the Doppler effect could lead, in principle, in the foreseeable future, to 
the measurement of the quadrupole moment of the giant planets of the Solar System. 

Keywords: Doppler Effect; General Relativiy. 



1. Introduction 

Einstein's theory of gravity, General Relativity (GR), has passed all observational 
tests with excellent results, at least at the scale of the Solar System and in many 
binary compact systems; however as it is well known, some problems arise with 
the recent cosmological observations. In fact, the evidence of the acceleration of 
the Universe is supported by experimental data deriving from different tests: i.e., 
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from Type-la Supernovae, from CMWB and from the large scale structure of the 
universe i' 2 * 3 ' 4 * 5 . However, GR is not able to provide a theoretical explanation of 
these experimental facts unless some exotic and invisible energy is admitted to exist 
in the universe (Dark Energy). Also for these reasons, the reliability of the theories of 
gravity alternative to GR is of great interest today. Of course, all these theories must 
agree with the known tests of gravity performed in the weak field and slow motion 
limit, in the Solar System: the standard tool used for formulating (metric) gravity 
theories in this context is the Parameterized Post-Newtonian (PPN) formalism 6 . 
Different theories are characterized by different values of the coefficients appearing 
in front of the post-Newtonian metric potentials, and the formalism determines 
the values of these parameters for each theory under consideration. Consequently, 
the formalism allows for a direct comparison of the competing theories with each 
other, and with the results of the experiments and observations. In particular, Solar 
System tests limit the range of the values that the PPN parameters may assume, 
thus excluding some theories of gravity and quantitatively limiting the deviations 
from GR 7 . 

The study of the Doppler effects in the space-time around a gravitating body was 
carried out by many authors in the past (see, for instance, 8 > 9 40,ii42 anc [ references 
therein), both for pedagogical purposes and for studying the feasibility of actual 
experiments aimed at detecting the effects themselves. In particular, in 11 , Kopeikin 
and Schafcr carefully studied the covariant theory of propagation of light in very 
general gravitational fields, and the various phenomena related, among which, the 
frequency shift. Indeed, as it is clearly explained in 13 , the Doppler frequency shift, 
the deflection of light and the time delay of electromagnetic signals are related to 
each other, in the sense that they might be thought of as "different manifestation 
of the same aspect of the [...] gravitational field". However, their measurements 
require different experimental techniques and, moreover, the interpretation of the 
measured quantities is quite different. In particular, the Doppler frequency shift is a 
rclativistic invariant, which, roughly speaking, corresponds to the projection of the 
photons' four-momentum along the world-line of the observer, and, consequently, it 
can be defined without reference to any particular coordinate frame. 

As for the measurement technique, in actual Doppler experiments a radio beam 
is transmitted from the Earth to the receiver located aboard a spacecraft; the radio 
beam is coherently transponded and then sent back to the Earth, where the received 
frequency is measured with great accuracy, usually by means of hydrogen masers. 
The comparison of the transmitted and received frequencies gives the measurement 
of the Doppler shift. 

In the past, there have been different proposals and different experimental tests 
of the relativistic frequency shift (see 10 and references therein), and, indeed, very 
recently, a test of GR was performed by studying the Doppler effect from the Cassini 
spacecraft. In particular, a measurement of the post-Newtonian parameter 7, which 
is equal to 1 in GR, was carried out: as a result of the experiment, 7 = 1 + (2.1 ± 
2.3) x 10~ 5 , thus sensibly limiting the deviations from GR 14 . 
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Here we apply the study of the Doppler frequency shift to some situations that 
may be of interest, both at the Solar System and astrophysical scale. In particular, 
we focus on Doppler effects originating from the bending of light rays in curved 
space-times (which will be measured, with great accuracy by the LATOR mission 
15 ), in order to see if these effects could be used to constrain the PPN parameters 
and, on the other hand, we argue that they can be used as probes to determine some 
astrophysical properties of interest of the celestial bodies. Numerical estimates of 
the magnitude of these effects are given, and they are compared to the current 
accuracies; furthermore, the characteristics of the effects are examined, in order to 
suggest possible detection techniques. 

The paper is organized as follows: in Section 2 we briefly outline the general 
framework used to study the Doppler frequency shift, according to previous works; 
in Section 3 we apply the general approach to some metrics of physical interest; in 
Section 4 we consider some prototypical situations and give numerical estimates for 
them; finally, the possibility of measuring the effects is discussed in Section 5. 



2. The Problem 

In this section we outline the theoretical framework which allows to study the 
frequency shifts of electromagnetic signals propagating in a gravitational field. In 
doing so, we follow the approach of Bertotti and Giampieri 9 , which is in agreement 
for the cases under consideration, with the very general approach developed by 
Kopcikin and Schafer 11 . 

Let us consider an emitter of monochromatic electromagnetic waves, which emits 
with proper frequency ve- The emitter is moving in the gravitational field of a 
massive source; let us call ue its four- velocity, evaluated at the (coordinate) time 
of emission of a signal when the emitter is at the position Se (the massive source 
is located at the origin of our reference frame). The problem is to evaluate the 
frequency vr measured by a receiver at xr, whose four- velocity is ur a . 

We may write the four- velocities in terms of the world- lines of the emitter (xe) 
and receiver (xr) 

dx E dx R 
ds E ds R 
Of course, in terms of the parameters se,sr along the world-lines: 

x e (se) = (x e (se),x e {se)) , xr(sr) = (x° r (sr),xr(sr)) . (2) 

Furthermore let £(A) = £o(A)+<!>£(A) be the world-line of the emitted photon, where 
£o is the trajectory that the photon would have if it were emitted in flat space-time 



a Greek indices run from to 3, Latin indices run from 1 to 3; the space-time metric has sig- 
nature (—1,1,1,1), and wc use units such that G=c=l; nonetheless, for the sake of clarity, we 
re-introduce physical units in the final formulae; boldface arrowed letters, refer to vectors in the 
three-dimensional space; boldface letters refer to four-vectors. 
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(i.e. a straight line) and <$£ is the perturbation induced by the gravitational field. 
This world-line is parameterized in such a way that the photon is emitted at Xe(se) 
when A = —Ie, and it is received at xr(sr) when A = Ir. 

Let us define the four-vector tangent to the photons' world-line 

^^ + ; ««A) 
dX dX dX 

where po is the flat space-time term and Sp is the perturbation induced by the 
gravitational field. If we set I = Ie + Ir the flat space-time term p can be written 
in the form 

P0 =(l,^^(l,n), (4) 

where n is a unit vector along the spatial path of the photons in flat space-time. 

The frequencies in the receiver and emitter reference frames are related by the 
following formula 16,17 

vr u R ■ k R 

— = r— , (5) 

v E u E ■ k E 

where k is the photons' wave vector. Indeed, in general, the frequency v p of a photon 
measured by an observer at a point P is found by projecting k onto the observer's 
four- velocity: 

2-KVp = kp ■ up. (6) 

Now, we can expand the components of the wave vector about their flat space-time 
values 10 : 

k = lo (p + Sp) = [-lo (1 + 5p°) , a (to + Sp)] . (7) 
Hence, eq. (5) can be written as 

VE \uE-P0)\l + y^)' 

where the fact that Po\e = Po\r has been exploited. From (8) we then obtain 

ln /M =ln (^Po\ +l n f i+ ^py (9) 

\v E ) \ue-PoJ \i + ^pe y> 

We see that the frequency shift (9) is made of two contributions. The first one 
depends on the velocities of the emitter and receiver and, also, on the gravitational 
field at the emission and reception points. In particular, it is possible to show 10 
that the ordinary Doppler effect (transversal and longitudinal) and the gravitational 
frequency shift come from this term. 

The second term depends, again, on the velocities and on the gravitational 
field, but, furthermore, it depends also on the gravitational perturbations of the 
world-lines of the photons. In the following, we focus on this term only. Also, we 
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assume that space-time is flat at infinity, which amounts to say that, in the limit 
for Ie,Ir —> oo, we may neglect the gravitational field at the emission and reception 
point. Consequently the first contribution accounts for the special relativistic effects 
(and, to lowest order, it is 0(v) h , where v is here the relative velocity of the emitter 
and the receiver), while the second term takes into account the gravitational effects 
along the world line. 

On neglecting the gravitational field both at the emission and reception points 
(indices are raised and lowered by means of the Minkowski tensor rj^), we may 
write the four- velocities of the emitter and the receiver 

u E /R = 1e/r (l^E/fl) , (10) 

where 

lE/R = ~7r = =5 = - 
V 1 ~ v e/r 

So, the contribution due to the perturbation of the photons' paths ("gravita- 
tional" contribution) reads 



In 



VR 



In 1 



dp — v ■ dp 
(1 — v ■ n) 



(12) 



where [z] R means z\r — z\e for an arbitrary function z. If we approximate the 
logarithm on the right hand side of (12), we may write 

ln f— ) =[(5p°-v-5p)(l + v-n)]*. (13) 

In order to calculate (13), we must evaluate the perturbations of the photons' 
momentum, at the emission and reception points, which can be done by integrating 
the geodesic equation. We reproduce here, without details, the results of 9 . 

The perturbations in the photons' four-momentum can be evaluated by solving 
the geodesic equations: 



d8p» 



K (p%+6p a )(pP + 5pe)=O. (14) 



We assume that the gravitational field is everywhere weak, so that it can be written 
in the form g M „ = r/^ + h^, where h^ v is treated as a perturbation of the Minkowski 
background -q^ v . Consequently, on using a linear approximation and supposing that 
the perturbations originate only in a small region of space around the source of 
gravitational field (this is the case if the impact parameter b is much smaller than 



b Since in our units c = 1, 0(v n ) means 0(v n /c n ), in other words we use approximations with 
respect to the small parameter v/c, where v is the scalar velocity. 
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E 



Z 



Fig. 1. Because of the presence of the massive source, the line linking the emission and 
reception points E, R bends, and does not follow a straight line. 

the distances of the emitter and receiver: b <C Ie, Ir), we can write 9 : 



The integrals (15), (16) must be evaluated along the unperturbed photons' path, 
i.e. along the straight line connecting the emission point and the reception point. 

Eq. (13), together with eqs. (15) and (16) are all we need to evaluate, to the 
desired order, the gravitational contribution to the Doppler frequency shift. 

3. Applications 

After having determined the perturbations of the photons' momentum, we may 
apply eq. (13) for evaluating the gravitational frequency shift. In the following, 
we study the frequency shift induced by a spherical distribution of mass, then we 
take into account a quadrupole contribution to the gravitational field and study the 
corresponding frequency shift and, finally, we study the effect of the rotation of the 
source of the gravitational field. 

3.1. The weak field around a spherical distribution of mass 

Let us consider a spherical mass distribution, and let us suppose that its gravita- 
tional field is weak. The line element can be written in the PPN form (up to 0(v 2 )) 




(15) 



(16) 
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18 

ds 2 = - (l - dt 2 + (dx 2 + dy 2 + dz 2 ) , (17) 

where r = \/ x 2 + y 2 + z 2 , and 7 is a PPN parameter. Let us suppose that 
photons propagate in the plane z = 0, and that their unperturbed path is parallel 
to the x axis, with impact parameter b c (see Figure 1). In other words, we have 

Po = (l,n) nil*. (18) 

Considering that we have no change of the time component of the photons' 
momentum, since the metric does not depend on time (" conservation of energy" ) , 
cq. (13), up to the lowest order in v becomes 

In = l(-v ■ 5p)f E = -v ■ 5p\ R + v ■ Sp\ E + 0(« 4 ). (19) 

Then, according to eqs. (15) and (16) 

Sp x \ E = S P x \ R = 0, 5p z \ E = 5p z \ R = 0, (20) 
and the only non null perturbations in the photons' world-line are given by 

6p v \e = - 1 4 T 2h xx ydx, (21) 



21 . 
I f°° 

5pV \ R = ii j s - h ( 22 ) 

On integrating (21) and (22) we obtain 

Sp^ E = 2 -^J^ l f Sp^ = - 2 -^^ l f. (23) 

As a consequence, the frequency shift due to a spherical distribution of mass, 
turns out to be 



In 



vr\ _2(1+ 7 )M 
ve 



A I 



(vr) v Ie + {v e ) v Ir 



+ 0(v 4 ). (24) 



3.1.1. On the physical meaning of the gravitational frequency shift 

What we have seen so far, allows us to give a simple interpretation of the gravita- 
tionally induced frequency shift. 

Let us consider an emitter moving with velocity v nearby a massive source (see 
Figure 2). Let us suppose that the emitter's velocity is orthogonal to the line con- 
necting it to the receiver. If we forget the gravitational field, we know that, in this 
case, there is no longitudinal Doppler effect (or first order Doppler effect, since it is 
0(v)), but there is only a transverse Doppler effect (which is 0(v 2 ), i.e. of second 



c In fact b is the ratio between two constants of the motion and coincides with the closest approach 
distance from the origin in the case of a fiat space-time. 



7, 2008 3:59 WSPC/INSTRUCTION FILE dopplerwfrev 



8 M.L. Ruggiero, A. Tartaglia, L. Iorio 




Fig. 2. The combined effect of the emitter's motion and of the gravitational bending of 
light rays produces a gravitational Doppler effect. 



order). This depends on the fact that the velocity of the emitter has no components 
along the line of sight. However, we know that one of the effects of the gravitational 
field is the bending of the light rays: this means that, roughly speaking, the light ray 
propagates along the curve 2 and not along the straight line 1. As a consequence, 
there is a component of the emitter's velocity along the actual direction of propa- 
gation of the signal, and its value can easily be calculated. In fact, the deflection 
angle 6 in GR (see, for instance, 16 ) is given by 

<») 



From Figure 2, it is easy to recognize that the component of the emitter's velocity 
along the propagation line is approximately given by 

4M 

v p ~ v20 = v ~^~- 

To this longitudinal velocity it corresponds a first order Doppler effect 

8v 



(26) 



VE 



or, reintroducing physical units, 



8v 



= v. D 



v 4GM 



c c 2 b 



(27) 



(28) 



which is in agreement with (24) if we set 7 = 1, as in GR, and let Ie <C Ir (see 
below). So, this frequency shift may be explained in terms of the bending of the 
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light rays due to the gravitational field, and it is a third order (0(v 3 )) Doppler 
effect. 



3.2. Quadrupole contribution in weak field approximation 

If the mass distribution is not perfectly spherical (still preserving an axial symme- 
try), the quadrupole contribution comes into play. In this case, the gravitational 
potential is written in the form 19 



#*0 



M 
r 



1- Ji 



R\ 2 3 cos 2 (9- 1 



(29) 



jfgz, contains 



where R is the average radius of the mass distribution, and J2 : 
the components of the quadrupole tensor Qij. 

In particular, the quadrupole contribution to the gravitational potential is 

M fR\ 2 3 cos 2 9 - 1 
<K*)Q = - J *\-) 2 • ( 30 ) 

Let us work out the quadrupole contribution to the gravitational Doppler effect. 
The line element is written in the form 



ds 1 



(1 + 2(f>) dt 2 + (1 - 2 7 0) (dx 2 + dy 2 + dz 2 ) 



(31) 



As before, we consider photons propagating in the equatorial plane (z = or 
9 = it/2) of the source, parallel to the x axis. This amounts to writing the metric 
perturbation h xx in the form 

2 

(h xx ) Q ■ (32) 



_ 2 7 M 7 M R\ 

flxx ~ ~ 2 ~^~ \~ ' -y h xx)M 



The first term, giving the monopole contribution to the frequency shift, has been 
calculated before; now we focus on the second term, expressing the quadrupole 
effect. 

Proceeding as before, we obtain that the only non null perturbations are given 

by 



Sp , lE . Wl + iiZ (|) is, V| . 



2J 2 (1+ 7 )M (R\ [e 
b J I 



Then, from 



In (j^j = {(-V ■ 60)]* - -v ■ Sp\ R + v ■ 5p\ E + 0(i; 4 ), 



(33) 



(34) 



taking into account the perturbations (33), we obtain for the quadrupole contribu- 
tion 

2 

+0(v 4 ). (35) 



R 



In = 2J 2(1+ 7 )M 



(vr) v Ie + (v E ) y l 



R 



I 
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3.3. The weak field around a rotating distribution of mass 

In the case of a rotating source and weak field approximation the line element, to 
the lowest meaningful order in PPN form, is 9 ' 7 

d S 2 = - (l - ^) ^ 2 +(l + 2J P) + V + ^ 2 )+ (47+ 2 ' 3 +Ql) (x x 3)^,(36) 

where S is the angular momentum and ot\ is another PPN parameter. As it is 
well known, the angular momentum contribution to the gravitational field, in this 
weak field context, is usually referred to as gravito -magnetic field 20 . 

Since the angular momentum contribution is of the order 0(v 3 ), we do expect 
that the gravito-magnetic contribution to the gravitational doppler shift is 0(v 4 ). 
In particular, it is possible to show 9 that this contribution is expressed by 




where ^5p means that we have to evaluate the perturbations up to 0(v 3 ). 

We consider, as before, photons propagating in the equatorial plane (z = or 
9 = 7r/2) of the source, parallel to the x axis; furthermore, we assume that the 
angular momentum of the source is parallel to the z axis: consequently there are no 
perturbations of the momentum along the x axis, and the non null perturbations 
are given by (see 9 ) 

( Vl^^(2(7+D + f), (38) 
( Vl«— ^(2(7+D + f). (39) 



So the gravito-magnetic contribution to the frequency shift turns out to be 



In 



(") 



S / . „, ai\ 

^( 2 ^ +1 ) + t) 



{vr) v Ie + {v e ) v Ir 



+ 0(v b ). 



(40) 



4. Examples and Numerical Estimates 

In order to evaluate the magnitude of the frequency shifts studied above, we consider 
three cases corresponding to the different positions of the emitter and receiver with 
respect to the source of the gravitational field: 

(1) Ie < Ir 

(2) Ir^Ie 

(3) l E - Ir 

The cases above correspond to different actual situations, that may prove useful 
for detecting the effect. Namely, case 1 can occur, say, in the Solar System, when an 
emitter is on board a spacecraft orbiting the Sun or one planet, and the receiver is 
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Source 


S (Kg m 2 s" 1 ) 


J 2 


Sun 


190.0 x 10 39 


2 x 10" 7 


Earth 


5.86 x 10 33 


1.08 x 10" 3 


Jupiter 


4.33 x 10 38 


1.49 x 10~ 2 


Saturn 


6.63 x 10 37 


1.62 x 10~ 2 



far away, for instance, on the Earth. Case 2 could correspond to a source of electro- 
magnetic signals, which is very distant from the massive source of the gravitational 
field. For instance, one may think of a distant star, such as a pulsar, emitting elec- 
tromagnetic signals, whose paths pass close to the Sun, and then are received on 
the Earth. Finally, case 3 may occur when both the emitter and the receiver are 
moving on similar orbits around a massive body. In the following, we shall evaluate 
the magnitude of the frequency shifts for some prototypical situations. 



Before going on, we notice that if we assume that 

vr = ve + Sv, (41) 

we may write 

ln (v_R\ =}n ( ^ + Su \/ + Su\Su 

\VE J V y E J V V E J VE 

since the expected frequency shift is small. 



4.1. Case 1., l E < I 



n 



As we said above, this situation may occur when the emitter is orbiting a massive 
source, and the receiver is far away from it. To fix the ideas, we may think of an 
emitter on a planet around the Sun, or on board a spacecraft orbiting the Sun or 
some other planets of the Solar System. Provided that 



Ie *C Ir, 



(43) 



we have Ir ~ /. As a consequence, the orders of magnitude of the various contri- 
butions that we have studied above may be evaluated by considering the following 
expressions, where, here and henceforth, for the sake of simplicity, the PPN param- 
eters are set equal to their GR values (i.e. 7 = 1, a\ =0): 



5v 



ve 
Sv 



ve 



M 



Q 



5v_ 

VE 



GM 



GM v E 
c 2 b c 

J 2 GM 
c 2 b 

GS v E 
c 3 b 2 c 



VE 

c 



(44) 
(45) 
(46) 
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Source 




iil 


M 




5v_ 

VE 


Q 




5v_ 

VE 


GM 


Sun 
Earth 
Jupiter 
Saturn 


3.06 x 10~ 9 
1.86 x 10~ 14 
2.28 x 10~ 12 
6.28 x 10~ 13 


6.12 x 10~ 16 
2.02 x 10" 17 
4.14 x 10~ 14 
1.02 x 10" 14 


1.41 x 10" 15 
3.03 x 10~ 20 
9.50 x 10~ 17 
1.32 x 10~ 17 



where ve, considering Keplerian motion, can be approximately written as 

V E - 



GM 



Rq 



(47) 



where Rq is the order of magnitude of the semi(major)-axis of the orbit. Conse- 
quently, we may write 

5v 



ve 
5v 



M 



VE 



GM GM 
^bi^Ro' 

J 2 GM (R 
c 2 b \~b 



GM 



VE 



GS GM 



GM 



More explicitly 



Sv 



VE 



~ 3,06 x 10" 



c 3 & 2 V c 2 Ro ' 



M_\ 3/2 (Re 



M 



Sv 



ve 
5v 



M 



MqJ 

3/2 



Mr. 



VE 



GM 



1,41 x 10" io I — 

bp, 



S 



3,06xlO- y J 2 [i- ^ [-- ^ 



R, 



R 



Rq 

2 



1/2 



Rq 



Rq 
Rq 



Rq 
b 



M 
Mq~ 



1/2 



Rq 
Rq 



1/2 



1/2 



(48) 
(49) 
(50) 

(51) 
(52) 
(53) 



The orders of magnitude of the frequency shifts (51), (52), (53) are evaluated in 
Table 2, where we have considered the Sun and some planets of the Solar System as 
sources of the gravitational field, and we have chosen orbits such that Ro ~ b ~ R. 



4.2. Case 2., l R < l E 

Another interesting situation corresponds to a source of electromagnetic signals, 
which is indeed very distant from the massive source of gravitational field. For 
instance, one may think of a distant star, such as a pulsar, emitting electromagnetic 
signals, whose paths pass close to the Sun, and then are received on the Earth. In 
this case, we have 

Ir < Ie, (54) 
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so that Ie — I- As a consequence, the velocity which plays a role is now the one 
of the receiver, i.e. the Earth-based observer. To fix the ideas, for a pulsar whose 
beam passes near the Sun, vr is nothing but the apparent velocity of the Sun in 
the sky, i.e. vr ~ 10 4 m/s. Then 

Sv GM v R 



5v 



M 



8v 



On evaluating for M 



VE 

M , S 
8v 



GM 
= Sq, 



c 2 b c 

J 2 GM 
~c 2 b~ 

GS v R 
c 3 b 2 c 



VR 

c 



VE 

5v 



we have 
~ 7.06 x 10" 11 



A I 



b 



VE 



7.06 x 10~ n J 2 



Q 



b 



5v 



VE 



GM 



(55) 
(56) 
(57) 

(58) 
(59) 
(60) 



It is clear that the most favourable position is near to the opposition, where the 
impact parameter is as close as possible to the radius of the Sun. Near to the 
opposition, we may pose 

db 



m = 



dt 



(61) 



As a consequence, the estimates of the frequency shifts can now be written as 



5v 



VE 

5v 



M 



GM 
~^b 



db 



cdt 



ve 



J 2 GM 





2 


db 


(C 


) 


cdt 



8v 



VE 



GM 



GS 



db 



cdt 



(62) 
(63) 
(64) 



We point out that the frequency shifts (62) and (64) were already obtained in a 
previous work 29 , using a different approach based on the calculation of the gravi- 
tational induced time delay ("Shapiro time delay"); furthermore, (62) corresponds 
to the recent measurements performed by Bertotti and collaborators 14 of the PPN 
7 parameter by means of radio ranging to the Cassini spacecraft. 



4.3. Case 3., l E 



R 



This situation may occur, for instance, when both the emitter and the receiver are 
orbiting the source of the gravitational field, and their orbits have similar dimen- 
sions. One might think of two satellites communicating with each other and orbiting 
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the Earth. To fix the ideas, one could consider two GPS satellites, even though in 
their present configuration they can communicate with the Earth-based stations 
only. However, the satellites of the forthcoming GALILEO positioning system, at 
least in their second generation, could be able to communicate with one another, 
so it is useful to evaluate the magnitude of the effect for this possible physical sit- 
uation. In order to evaluate the magnitude of the effect, we consider two satellites 
orbiting the Earth, on the same circular orbit with radius R c ~ 26.6 x 10 6 m; wc 
also assume that the orbit lays in the equatorial plane of the Earth. Furthermore, 
we calculate the frequency shift in the most favourable position of the two satellites, 
i.e. when b ~ R®. As a consequence, we may write 
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Numerically evaluating the formulas, we obtain 
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(65) 
(66) 
(67) 

(68) 
(69) 
(70) 



GM 



As a general comment, we may say that, in all cases analyzed above, the gravito- 
magnetic contribution to the frequency shift is much lower than the other ones, so 
it would hardly be detected. Also, the value of the J2 (Table 1) for the bodies of the 
Solar System, makes the quadrupole contribution much smaller than the monopole 
one. 



5. Discussion and Conclusions 

We have shown that the relative frequency shifts originating from the gravitational 
bending of light rays range from 10~ 9 to 10~ 22 . The measurements of these fre- 
quency shifts, could allow, in principle, to estimate the PPN parameters that we 
have introduced in our formulas, or to measure the gravitational field of the celestial 
bodies. Actually, as we pointed out before, a recent estimate of the 7 parameter, 
was performed, by means of the measurement of the frequency shift of radio photons 
to and from the Cassini spacecraft as they passed near the Sun 14 . We remember 
that, in GR, 7 = 1, ct\ = 0, and the most recent values of the PPN parameters arc 
reported in 7 . 
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Of course, in order to evaluate the reliability of actual measurements of these 
effects, a careful and detailed analysis of the error budget is required. In this work, 
however, we are just interested in pointing out the possibility of detection of some 
of them, on the bases of the knowledge of the current accuracy in frequency shift 
measurements. 

We see that, in order to measure the quadrupole contribution, an accuracy of 
10~ 14 at least is required, while the gravito-magnetic contribution requires an ac- 
curacy of 10~ 15 . Nowadays, it is known that an atomic standard stable to 1 part in 
10 16 is available 30 , so that such a fractional frequency stability is not unreasonable 
for future experiments with atomic clocks. Indeed, the giant planets, Jupiter and 
Saturn, are, at least in principle, candidates for such measurements: in fact, the 
estimates of their quadrupole moment seem to be in the actual or foreseeable accu- 
racy of the instruments. Furthermore, they do not have an electromagnetic activity 
as intense as the one of the Sun: so, again, at least in principle, it would not bias 
the radio communications. 

We must remember that in an actual measurement, besides the gravitationally 
induced terms, also the ordinary Doppler effect (0(vj), the second order Doppler 
effect (0(v 2 ) and the gravitational redshift (0(v 2 )) contribute to the Doppler sig- 
nal (in particular, the latter two effects were tested using a space-borne hydrogen 
maser carried aboard a rocket 31 ). In order to measure the frequency shifts induced 
by the gravitational bending of light rays, these "ordinary" Doppler contributions 
must be subtracted: to this end, a very accurate knowledge of the orbits of the 
spacecrafts or planets is required. On the other hand, the different dependence on 
the geometric parameters (such as b) of the various contributions, may suggest a 
way of discriminating them; in the case of the gravito-magnetic contribution to the 
frequency shift, the antisymmetry of the effect with respect to the rotation axis of 
the source may help to disentangle the extremely weak signal from the remaining 
bigger shifts due to other causes 29 . 

In conclusion, we have shown that the Doppler effects induced by the gravita- 
tional bending of light rays, though small, deserve further analysis. If measured, 
these effects may allow us to put limits on the values of the PPN parameters and, 
on the other hand, they can be used as probes to understand the nature and explore 
the physical properties of the celestial bodies. 
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